Generalized correlation functions for conductance fluctuations and the mesoscopic spin Hall effect 
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We study the spin-Hall conductance fluctuations in ballistic mesoscopic systems. We obtain universal ex- 
pressions for the spin and charge current fluctuations, cast in terms of current-current autocorrelation functions. 
We show that the latter are conveniently parametrized as deformed Lorentzian shape lines, functions of an ex- 
ternal applied magnetic field and the Fermi energy. We find that the charge current fluctuations show quite 
unique statistical features at the symplectic-unitary crossover regime. Our findings are based on an evaluation 
of the generalized transmission coefficients correlation functions within the stub model and are amenable to 
experimental test. 
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I. INTRODUCTION 

The discovery of the spin Hall effect (SHE)A^ in both 
metal and semiconductor structures has opened an important 
new possibility to control the effects of non-equilibrium spin 
accumulation. 5 The basic idea underlying the SHE is to gen- 
erate spin currents transverse to the longitudinal transport of 
charge by creating an imbalance between the spin up and spin 
down states^ 

The detection of spin-Hall conductance fluctuations is a 
major goal of semiconductor spintronics. 8 It is, however, a 
hard endeavor. The main reason is the difficulty to efficiently 
connect ferromagnets leads to two-dimensional semiconduc- 
tor structures,^ 1 For ballistic systems coupled metallic leads, it 
is in principle possible to detect the signal when scattering by 
impurities induce a separation of the spin states. 

Through this mechanism, universal spin-Hall conductance 
fluctuations (USCF) can lead to accumulation of spin at the 
electron reservoirs. The USCF appear in the transverse cur- 
rent measured in multi-terminal devices in the presence of a 
sufficiently large magnetic field. 6 Signals of the spin accumu- 
lation can be inferred, for instance, from time-dependent fluc- 
tuations of the spectral currents (noise power)iS or from the 
analysis of universal conductance fluctuations Mr— Spin-Hall 
conductance fluctuations have been theoretically studied for 
mesoscopic systems in the diffusive-!! as well as in the bal- 
listic regime. 12 In the absence of both spin rotation symmetry 
and magnetic field, these studies predict universal spin-Hall 
conductance fluctuations with a root mean square amplitude 
of about 0.18(e/47r). 

So far, a direct detection of spin-Hall currents by analyz- 
ing transverse current fluctuations has not been made. In this 
paper, we propose an alternative way to infer spin-Hall con- 
ductance fluctuations, based on the universal relation between 
spin and charge current fluctuations in chaotic quantum dots. 
We find that the change and spin current-current correlation 
functions show a quite unique dependence on the ratio of open 
modes between transversal and longitudinal terminals. This 



dependence allows one to infer the magnitude of the spin cur- 
rent. From the technical point of view, we adapt the diagram- 
matic technique developed to describe the electronic transport 
in two-terminal chaotic quantum dots in the presence of a 
spin-orbit interaction, at the symplectic-unitary corssoverjA 5 *^ 
to the case of multi-terminal spin resolved currents. 

The paper is organized as follows. In Sec. [II] we review 
the Landauer-Biittiker approach used to calculate the multi- 
terminal charge and spin currents. Next, in Sec. [Ill] we present 
the diagrammatic theory we employ to calculate the universal 
spin Hall current fluctuations. The phenomenological impli- 
cations of our findings are discussed in Sec. [IV] Finally, in 
Sec.[V] we present our conclusions. 



II. THEORETICAL FRAMEWORK 

In this section we review the scattering matrix formalism 
that describes the spin Hall effect in ballistic conductors. We 
follow the approach put forward in Refs. [l2lll7lll8l We find 
helpful for the reader to have, in a nutshell, the expressions for 
the charge and spin currents with their explicit dependence on 
the electron energy and the presence of an external magnetic 
field. The latter is of particular importance in our study: The 
magnetic field breaks time-reversal symmetry and drives the 
symplectic-unitary crossover. The description of the univer- 
sal spin current fluctuations at the crossover is one of the key 
results of this paper, as discussed in SecHTH 

We consider multi-terminal two-dimensional systems 
where the electrons flow under the influence of a spin-orbit 
interaction of the Rashba or/and Dresselhaus type. The de- 
vice is connected by ideal point contacts to N independent 
electronic reservoirs, denoted by i = 1, • • • ,Af. The elec- 
trodes are subjected to voltages denoted by V*. We use the 
Landauer-Biittiker approach to write the a direction spin re- 
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solved current through the ith terminal as^ 

^ = jT,Y, \ S ^n.,'\ 2 (V t - Vj) (1) 

where a and a' are the spin projections in the a = x, y or z 
direction and S is the quaternionic scattering matrix of order 
2Nt x 2Nt that describes the transport of electrons through 
the system. The total number of open orbital scattering chan- 
nels is Nt = Y^iLi where JV* is the number of open chan- 
nels in ith lead point contact. 

The electric current at the ith terminal is i> ' = if^ + 1^, 
for any direction a = x, y, z of the electron spin projection. 
Similarly, the a-axis component of the spin current 1^ is de- 
fined as the difference between the two spin projections along 
the a-axis, namely, 1^ = 1"^ — 1"^, 

Let us consider a set up with TV = 4 terminals. An applied 
bias voltage between electrodes i = 1 and 2 gives rise to a lon- 
gitudinal electronic current / and due to the spin Hall effect 
to spin currents at the transversal contacts.— Charge conser- 
vation imposes that / = if* 1 = —1%. Moreover, due to the 

absence of a transversal voltage bias, 1^ = for i = 3, 4. 
Using these constraints, it was showni^ that the transversal 
spin currents are given by 

4 a) = \ (li a) - 7t ] ) - T^V 3 - 1«>V A (2) 

where i — 3, 4. Likewise, the longitudinal charge current 
reads 

jf ) = i (2 Nl - rA 0) + 2N 2 r 2 ( 2 0) + r x ( 2 0) + r 2 { ?) 

-l^? -7^)Vs-i(7i 0) -7#>)fc, (3) 

where i = 1,2. For notational convenience one introduces 
the dimensionless currents J = h/e 2 (I /V). The effective 
voltages Vt (in units of V), are given by rather lengthy expres- 
sions of the generalized transmission coefficients T~j , which 
can be found in Ref. Il2j Finally, the generalized transmission 
coefficients read 

7^\e, E', B, B') = tr [(L ® a a ) &(E, B)ljS(E', B')] , 

(4) 

where E and B stand for the electron energy and the mag- 
nitude of an external magnetic field. The trace is taken over 
the scattering channels that belong to the ith and jth point 
contact. The matrix L stands for a projector of the scatter- 
ing amplitudes into ith point contact channels. The matrices 
cr Q , with a £ {x, y, z}, are the Pauli matrices, while a is the 
2x2 identity matrix characterizing an unpolarized (charge) 
transport. 

III. UNIVERSAL SPIN-HALL CONDUCTANCE 
FLUCTUATIONS 

We assume that the electron dynamics in the quantum dot 
is ballistic and ergodic^Si 2 ! and model the system statistical 



properties using the random matrix theory. In this section, 
we describe the procedure to obtain universal spin and charge 
current-current autocorrelation functions given by Eqs. (0 
and ([3). We compare our analytical expressions the with re- 
sults from numerical simulations. 

Let us assume that spin-orbit interaction is sufficiently 
strong to fully break the system spin rotational symmetry. In 
other words, the spin-orbit scattering time is much smaller 
than the electron dwell time in the system, that is, r s . . -C 
Tdwcii- Hence, in the absence of an external magnetic field, 
the scattering matrix has symplectic symmetry. By increasing 
B, the system is driven through a crossover from the symplec- 
tic to the unitary symmetry. 

The resonance S'-matrix can be parameterized as 25 

S(E,B) =TU[l-Q^R(E,B)QU]~ 1 T\ (5) 

where U is a matrix of order 2M x 2AI of quaternionic form, 2 - 
M stands for the number of resonances of the quantum dot. 
We take M ^> Nt- The matrices Q and T describe projector 
operators of order (2M - 2N T ) x 2M and 2N T x 2M, re- 
spectively. Their matrix elements read Qjj = 5i+2N T ,j and 
Ti,j = Sij. The matrix R of order (2M-2N T )x(2M-2N T ) 
reads 

R(e, x) = cxp (^<r a + ^Xa ^ , (6) 

where X is an anti-hermitian Gaussian distributed random 
matrix, while e = (Ta wc \\/h)E and x 2 = t^qw/tb are di- 
mensionless parameters representing E and B. 

Both the dwell time Tdwcii and the magnetic scattering time 
tb are a system specific quantities. Tdwdi is usually ex- 
pressed in terms of the decay width, namely T = h/T^ wc \\ = 
NtA/2tt, where A is the system mean level spacing. In 
turn, Tg 1 is the rate by which the electron trajectory accumu- 
lates magnetic flux in the quantum dot. For chaotic systems, 
Tg 1 = k(AB) / 4>o, where <po is the unit flux quantum, A is 
the quantum dot lithographic area, and k is a diffusion coeffi- 
cient that depends on the quantum dot geometry r— 

The S-matrix given by Eq. (0) can be formally expanded in 
powers of U, namely, 

oo 

S(e,x) = J2 TU [Q^R{^x)QU] m T\ 
m=0 

and used to calculate the ensemble average of the transmission 
coefficients {T^) for M 3> Nt Following this algebraic 
procedure, we write the sample transmission coefficient as a 
two-point function of the S'-matrix, 

oo 

Tf\e,e',x,x') = £ Tr{(L ® a a ) T [U^Q^R\e, x)Q] m 

m,n— 

x U^T1 3 TU [Q f -R(e', x')QU] "r 1 }. 

(7) 

For chaotic quantum dots, as standard,— we assume the ma- 
trix elements of U as Gaussian variables, with zero mean and 
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variance 1/M. This allows one to express the calculation of Tr (lj <g> a a ) = 2Ni5 a o and Tr (lj) = 2Nj to write 



moments and cumulants of 7y by an integration over the 
unitary group leading to a diagrammatic expansion in powers 
of 1 in terms of diffusons (ladder) and cooperons (max- with 
imally crossed) diagrams. The method is described for the 
Dyson ensembles in Ref. 25. This approach has been ex- 
tended to treat the crossover between symmetry classe o 15 ' 27 
and is shown to render the same results as the quantum circuit 
theory M 

The diagrammatic technique 2 - allows one to calculate mo- 
ments up to arbitrary order of the transmission coefficients. 
It consists in grouping and integrating over the Haar measure 
all the independent powers of U in Eq. (01. Figure Q] shows 
the diagrams that represent the leading order contributions to 
all possible contractions of the U matrices. In the sum of Eq. 
01, we verify that, after taking the average over U, only the 
powers with m = n contribute to the average transmission 
coefficients. The white and black dots in Fig. Q] stand for the 
indices of the matrix U, with elements E/Jy in the channels 



space, while the Greek symbols represent the Pauli matrices 
indices a p . a in the spin space. 



a' p' 



a 1 a' p p 



a p 



FIG. 1 : Diagrams representing the leading contributions to the aver- 
age transmission coefficient (7y ) of Eq. ((7}. Vertical lines repre- 
sent contractions of the U matrix elements. The white and black dots 
in stand for the channel indices of the matrix U and the Greek sym- 
bols represent the Pauli matrices indices a Pt<7 in the spin space. The 
first diagram (top) is known as diffuson and accounts for the semi- 
classical diffusive processes. The two others diagrams are known as 
cooperons and give the mean quantum correction, namely, the weak 
(anti)localization corrections. 



The diagram at the top row of Fig. Q]is called diffuson, in 
analogy to the ladder diagram that appears in diagrammatic 
expansion to calculate the conductivity in disordered diffusive 
mesoscopic systems^ The diffuson contribution to the aver- 
age transmission coefficient reads 



(7#° (e, e', x, x')) (d) =\ £ [Tr (1 4 ® a a ) OTr (1 

per 

-Tr(l 4 ®(i")2? 2 Tr(l J ) 



pa\a p 



(8) 



where, for clarity, we make explicit the spin degree free- 
dom of the symplectic structure. We use the properties 



V- 1 = 2Ma° ® cr° - Tr (R ® R' 



Tr (R <g> fit) = a ® ct° 
x [2M - 2N T - 2i(e 



(9) 



(10) 



x>) 



/\21 



where R' is a shorthand notation for R(e',x'). The tensor 
products follow the "backward algebra" ) ' 6 namely, (a 1 <£> 
(j j )(a k ® a 1 ) = (a i a k ) <g> (a l a^). 

The diffuson contribution to the generalized transmission 
coefficient is obtained by evaluating Eq. ([8} using the expres- 
sions (O and ( [Tol l. It reads 



{7-f( e , £ ',x,x')} 



(d) 



1 



TT- , (ID 



xf/2}. 



where N D = N T [l - i (e' - e) + {x 1 

We are now ready to evaluate the two maximally crossed 
diagrams at the bottom row of Fig. Q] They are known as 
cooperons 2 ^ and represent the main quantum interference cor- 
rection to the conductance in chaotic systems, responsible for 
the weak localization peak. 

Following the procedure described above, we obtain the 
cooperon contribution for the generalized average average 
transmission, namely, 

(T^\e,e',x,x')}^ = ^ \M F ^ {TJttT)F. 



(2M> 



rTrfi^ (Tfrr7i)Tr{F s 



(12) 



The operator T = cr° ® <J V is related with the time-reversed 
of path in the cooperon channel of the dual space. We also 
define the matrices 



F; 



(«) 



h 



cr Q +Tr(l i (g)(T")P(i?' t 
■Ti(lj)V(R* ®R), 



R), 



where 



fuu = [2M<r ® a - Tr(i? $ R*)] 1 
f TT = (2Ma° ® cr ) Tr(R ® R'*)fuu- 



(13) 



(14) 



The superscript * denotes the quaternion complex conjuga- 
tion. Using the quaternionic conjugation rules and taking the 
limit M — > oo, we obtain 



/^=2iV c a ® ( T . 

f TT = (2Ma Q ® (T )[(2M - 2N C )<J (8 <T ]fuu, 



(15) 



where N c = N T [1 -i(e' -e) + (x' + x) 2 /2}. 

Summing up the diffuson and the cooperon contributions to 
the generalized transmission coefficients we obtain 



(Vi ,x,x )} = S a0 <{ J 



N c 

NjiNp+ije-e')] 



(16) 
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For a = 0, Eq. ( [Tol l reproduces the average electron transmis- 
sion reported in the literature^. As expected in this case, the 
average spin transmission coefficients are zero. 

Let us now analyze the transmission coefficient fluc- 
tuations. We use the same diagramatic procedure de- 
scribed above for all the 32 diagrams characteristic of the 
usual covariance calculations.— To address relevant phys- 



ical situations, it is sufficient to consider transmission 
coefficients with single energy and magnetic field argu- 
ments, that is, T^{e,x) = 7^- (e, e, x, x). Following 
the diagrammatic approach^ we calculate the covariance 

cov[7Tf (e,x),T^(e',x')} ee (7$*\e, x)T^(e', x')) - 
(T^\e,x)}(T^\e',x'))^d obtain 



T^\e,x),TlP{^x') 
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FIG. 2: (Color online) Generalized average transmission coefficient 
(7^°^ (e, x = 0)) versus energy e = E/Y for different spin and ter- 
minal indices. The analytical results (stub model) are represented by 
the solid and doted lines, while the results of numerical simulations 
(RMT Hamiltonian model) are represented by the symbols. The fig- 
ure inset describes the different 7"' 
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Support to our analytical findings is provided by numeri- 
cal simulations. For that purpose, we find convenient to em- 
ploy the Hamiltonian approach to the S-matrix^. The lat- 
ter is more amenable for numerical simulations than the S 
matrix parametrization of Eq. ([5]) and both are statistically 
equivalent^ 

The Hamiltonian parametrization of the S matrix reads 

S{E,X) = l-2niW^(E-H(B)+iwWW^)~ 1 W , (18) 

where E is the electron propagation energy and H(B) is the 
matrix of dimension 2M x 2M that describes the resonant 
states (M orbital states times the 2 spin projections). In gen- 
eral, H depends on one (or more) external parameters X. As 
discussed before, we are interested in the case where by in- 
creasing B one breaks time-reversal symmetry, driving the 
system from the symplectic to the unitary symmetry. In our 



numerically simulations, we consider B 3> B c , namely, the 
specific case of pure ensembles. Accordingly, H is taken as 
a member of the Gaussian unitary ensemble corresponding to 
the case of broken time-reversal symmetric case, usually de- 
noted by /3 = 2. The matrix W of dimension M x (2Nt) con- 
tains the channel-resonance coupling matrix elements. Since 
the H matrix is statistically invariant under unitary transfor- 
mations, the statistical properties of S depend only on the 
mean resonance spacing A, determined by H, and W^W . We 
assume a perfect coupling between channels and resonances, 
which corresponds to maximizing the average transmission 
following a procedure described in Ref . [3jl 

For simplicity, we take the case of N = iV,. The results 
presented in Figs. [2] and [3] correspond to the systems with 
N = 5 perfectly coupled modes and M = 400 resonant 
levels. Hence, the S'-matrix has 2Nt = 40 open channels. 
The ensemble averages are taken over N r = 10 5 realizations 
within an energy interval around the band center, comprising 
about M/A resonances. 

Figure [2] compares the average transmission (7y (e)) ob- 
tained from the numerical simulations with the analytical ex- 
pression ( fToT ) for a number of different cases. The agreement 

— 1/2 

is very good, with accuracy of the order of N r ' . The sim- 
ulations indicate that the average transmission in stationary in 
e = E/T, as it should. 

Figure [3] contrasts transmission coefficient covariances cal- 
culated using Eq.dTTb with those obtained from numerical 
simulations for a number of different cases. As before, the 
discrepancies are very small and stay within the statistical pre- 
cision N r 1 ^ 2 . The random matrix theory^ predicts an auto- 
correlation length r = iV^A/ (2ir) for a two-terminal geome- 
try. Our results for the correlation function extend the latter to 
four-terminal geometries with (or without) spin polarization. 

Let us now return to the problem of spin and charge cur- 
rent and effective potential. As mentioned, both are combi- 
nations of transmission coefficients. Fortunately, it is possi- 
ble to calculate average currents and current-current correla- 
tion functions in terms of the average and the transmission 
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FIG. 3: (Color online) Transmission coefficient covariance 
cov\T^"\e,x = 0)7^; (e', x = 0)] as a function of the energy 
difference Ae = (E — E') /T for different spin and terminal indices. 
The analytical results (stub model) are represented by the solid and 
doted lines, while the results of numerical simulations (RMT Hamil- 
tonian model) are represented by the symbols. The figure inset de- 
scribes the different cav\T{"\e)T^P {e')\ considered. 



coefficients correlation functions already calculated and con- 
firmed numerically. The effective voltages V3 and V4 show 
sample-to-sample fluctuations that depend both on the en ergy 
and magnetic field. On the other hand, as discussed in Ref.[l2l 
their ensemble averages depend only on the number of open 



channels, namely, {Vi(e,x}} = 1/2(N 1 - N 2 )/{N l + N 2 ), 
with i = 3,4, We also note that the ensemble average of the 
spin current is always zero, ( J> a ' (e, x)) = 0, with i = 3,4 
and a^0, independently of the energy and magnetic field. 

The USCF do not depend on the device geometry (nor on 
the positions of the terminals), but rather on the number of 
open channels at each terminal. Without loss of generality, 
let us analyze the spin current covariance for the case N\ = 
N2 = N and N3 = N4 = nN, a setting that is easily realized 
in experiments. Here, n is a real positive number that we call 
"channel factor"^ 

For the spin currents, for which a ^ 0, we obtain 



cov 



jf Q) ( e ,zU a) ( e V) 



1 n/{l + n) 2 
8 (1 + Sx 2 ) 2 + Se 2 



(19) 



where Se = e — e' and 5x = x — x' . It is worth noticing 
that, for n = 1 and Se = Sx = 0, Eq. ( fT9b perfectly repro- 
duces the recent reported result o 11 ! 12 for the universal fluctua- 
tions of the transverse spin conductance, namely, rms[G s n] = 
e/4n{cov[J i l a) (e,x),J^ ) (e,x)}} 1 / 2 w 0.18e/47r. 

Equation dl~9b shows that the spin current correlation func- 
tions do not depend on the cooperon channels, that give rise 
to terms containing Nc- Hence, these quantities do not de- 
pend on the magnetic field, represented by x, but rather on its 
variations, Sx. As a consequence, in the set up we consider, 
the spin current fluctuations are invariant in the symplectic- 
unitary crossover regime, a quite remarkable property. 

The charge current fluctuations, on the other hand, depend 
both on the cooperon and diffuson channels, leading to 



cov[4°\e,x),jP(e',x')} = ±-{ ^ + 2n 1 + 1 \ — — - 1 , (20) 

16 y (1 + n) 6 (1 + Sx 2 ) + Se 2 (1 + n) A [I + (2x + Sx) 2 ] 1 + Se 2 J 



where i = 1,2. The magnetic field, represented by x, drives 
the symplectic-unitary crossover. For x = 0, one recovers 
the symplectic limit, while the unitary one is attained when 
x ;§> 1. Note that in the absence of "transverse" leads, or 
n = 0, Eq. (l20l reproduces the two-terminal result found in 
Ref.Ql 

From Eq. ( f20b it follows that 



cov{[jf\e,x) 



4 0) (e, 



16nx(l + 2x 2 ) 
(l + n)(l+4a; 2 ) 5 



(21) 

which demonstrates that, except for the two-terminal case 
where n / 0, the charge currents are not even functions of 
the magnetic field x<^ 



IV. ALTERNATIVE STATISTICAL MEASURES 



Equations (fT9l and (l2Qb are the main results of this paper. 
Unfortunately, the statistical sampling required to confirm our 
predictions for the dimensionless currents is rather large, mak- 
ing the experimental requirements quite daunting. An easier 
accessible statistical measure has been recently proposed: 2 ^ 
The dimensionless current J^' fluctuates as e and x are var- 
ied. Let us call the external parameter z. Useful statistical 
information can be extracted from the number of maxima (or 
minima) of the j[ (z) in a given interval [z,z + Sz]. Us- 
ing a scale invariance and maximum entropy principle, we re- 
late the joint probability of j\ (z) and its derivatives to a 
general equation for the density of maxima, for spin and/or 
charge transport. The average densities of maxima, (ftz ) of 
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the fluctuating current j[ are given by 20 

»\ = 

v[t\e, X ),j\ a \e', X ')\ 



T 2 



'27rV T 2 

d 2 

d{5z) 2 
d 4 



Sz=a 



d(5z) 4 



Sz=0 



where 5z is 5e or Sx, 

It is convenient to write the charge current covariance as 
a deformed Lorentzian. For parametric variations of z, we 

setcov[4 0) (e,x),jP(e,x')} = a z (n,x)/[l + {5x) 2 ] h ^ x \ 
where a z (n, x) is a crossover function and h z (n, x) charac- 
terizes the Lorentzian shape deformation of the charge current 
correlation as a function of Sx. In terms of the factor h z , the 
average density of maxima reads 



(23) 



where z can be either e of x. 




FIG. 4: [Color online] Contour plot of h x (n, x) as a function of the 
magnetic field, represented by x, and the channel factor n. The color 
code is explained at the strip on the right panel. 

Using Eqs. d20l ) and d23l . we obtain an exact analyti- 
cal expression for h x (n,x). Its explicit form is not pre- 
sented here, since it is rather lengthy. Figure [4] illustrates 
its general features. For the unitary symmetry limit, it is 
well stablishe d 34 ' 35 that the electronic conductance correlation 
function shows a square Lorentzian behavior Accordingly, we 
find h sqL = lim-r^oo h x (n,x) = lim x _>o h x (n, x) = 2, for 
the pure circular unitary and symplectic ensembles, respec- 
tively. The symplectic-unitary crossover shows a much richer 
behavior. Figure [4] exhibits a remarkable crossover between 
sub-Lorentzian, for which h x < 1, with a minimum value of 
h x « 0.64, and super-Lorentzian, for which h x > 1 with a 
maximum value of h. x « 2.92. 
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FIG. 5: [Color online] Contour plot of h t (n, x) as a function of the 
magnetic field, represented by x, and the channel factor n. The color 
code is explained at the strip on the right panel. 



Parametric variations of e were first studied in nuclear scat- 
tering at low energies and known as Ericson fluctuations 36 . 
As it is well-known, their characteristic correlation func- 
tion versus energy has a Lorentzian shape. In the pres- 
ence of a perpendicular magnetic field and the channel fac- 
tor, we obtained a unitary-simplectic crossover of Lorentzian- 
type shapes, generalizing the correlation functions of Er- 
icson fluctuations. For parametric variations of e, we set 



cov[J < w (e,as),4 UJ (e' J x)] = a e {n,x)/[l + (8e) 2 ] h ^ x \ 
where a e (n,x) is a crossover parameter and h e (n,x) char- 
acterizes the deformation of the Lorentzian shape. As in the 
previous case, we also obtain a lengthy analytical expression 
for h e (n,x). Its main features are displayed in Fig. [5] As 
expected, = ]im x ^ too h e {n,x) = lim. E _ i .o h e (n, a;)=l, for 
pure circular unitary and symplectic ensembles, respectively. 
Figure [5] exhibits another remarkable crossover from a sub- 
Lorentzian, for which h e < 1, to a Lorentzian behavior. 

According to Eq. d23l . the density of maxima correspond- 
ing to pure ensembles, namely, x = or x 3> 1, is (px) = 
y/6(h^ L + 1) Ps 0.68 and (p { e a) ) = y/6(h L + 1) w 
0.5522iSi, for both spin and charge currents. We emphasize 
that for the case of the spin correlation function, Eq. $1% , 
h c = hi, = 1 and h x = /i sq L = 2 even in the crossover 
regime (any value of n and x). 

Let us now focus on the longitudinal (charge) correla- 
tion function, Eq. ( f20b . For a given value of the channel 
asymmetry factor, (p x ^ (n, x)) has a unique global maximum, 
(pi (n, x max )), and minimum, (pi'(n, x min )) The differ- 
ence, A(p x (n)) = (pi (n, x rnax )) - (p x °\n, x rnin )) in- 
creases with n until it saturates at n ~ 3. In the absence 
of spin leads, we find the difference A(p x (0)) ps 0.27. 
In the presence of spin leads, we get A (p x (0.5)) ~ 0.19, 
A ( Px (l)} « 0.17, A (^(2)) « 0.15, and A (p x (5)) « 0.14. 
Thus, in measurements made with a perpendicular magnetic 
field, for symmetric channels (n = 1), the spin terminals lead 
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to a reduction in the signal of about 37%, which becomes 
even larger with increasing n. Interestingly, the maximum 
and minimum of (p x ) correspond to magnetic field strengths 
Xmin ~ 0.20 and x max ps 0.47, for n S [0, 5], a rather narrow 
range of values which is accessible experimentally. 

In contrast to (p x {n)), the energy variation generates a den- 
sity of conductance peaks containing a single global mini- 
mum, (pe(n, x m i n )), andno global maximum. This minimum 
is located in a very narrow range of values of x m i n « 0.26 as 
a function of n. The minimum value of the density at these 
points is of the order of (p e ) m 0.52. 

V. SUMMARY AND CONCLUSION 

In this paper, we have investigated the spin-Hall conduc- 
tance fluctuations in a chaotic open quantum dot with spin- 
orbit interaction. Both the electronic and the spin-Hall con- 



ductance fluctuations are universal functions, with autocorre- 
lation functions that depend on the magnitude of the exter- 
nal magnetic field B and the channel asymmetry factor n. A 
clear intermediate case of symplectic-unitary transitional be- 
havior is found and can be tested experimentally. In particu- 
lar, the spin current can be measured by using the charge cur- 
rent density of maxima. The results of this Letter extend the 
understanding of mesoscopic fluctuations to spin- and charge 
currents in the symplectic-unitary crossover, characteristic of 
quantum dots subjected to an external magnetic field. 
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